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Abstract 

In this paper we consider FRW cosmology in modified non-local gravity. The stability analysis 
shows that there is only one stable critical point for the model and the universe undergoes a 
quintessence dominated era. 
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1. INTRODUCTION 



Cosmological observations such as Super- Nova la (SNIa), Wilkinson Microwave 
Anisotropy Probe (WMAP), Sloan Digital Sky Survey (SDSSV, Chandra X-ray Observa- 
tory disclose some cross-checked information of our universe jl|-j4|. Based on them, the 
universe is spatially flat, and consists of approximately 70% dark energy (DE) with nega- 
tive pressure, 30% dust matter (cold dark matters plus baryons), and negligible radiation. 
The combined analysis of SNIa, that is based upon the background expansion history of 
the universe around the redshift z < 0{1), galaxy clusters measurements and WMAP data, 
offers an evidence for the accelerated cosmic expansion jsljo]. The cosmological acceleration 
strongly indicates that the present day universe is dominated by smoothly distributed slowly 
varying DE component The ACDM model with an equation of state (EoS) parameter 
being —1 has been continuously favored from observations. From Cosmic Microwave Back- 
ground (CMB) and Baryon Acoustic Oscillations (BAO) jsljo], a new constraint on the EoS 
parameter are observed to be around the cosmological constant value —1 ±0.1 {g]- 10| and 
a possibility that the parameter is time dependent. 

The DE models with the possible phantom crossing can be broadly classified into two 
categories, (I) by adding one or more scalar field into the standard formalism of gravity, for 



example see ( ll|-tl6 and refs. therein) (II) by any kind of modification in the geometry 
of the gravity, see |l7|-j23|. However, a major issue in analyzing cosmological models stems 
from the fact that the field equations are nonlinear and thus limits the possibility of obtaining 
exact solutions and analyzing the behavior of the cosmological parameters. On the other 
hand, in recent years, an increasing realization of the importance of the asymptotic behavior 
of cosmological models in comparison with the observational data, emphasis study of the 
qualitative properties of the equations and of the long-term behavior of their solutions 



29| with different approach, we 



In this context, similar to the work by the authors in 
use perturbation method, especially linear stability and phase-plane analysis to study the 
stability of a modified non local cosmological model and investigate the possible quintessence 
scenarios in the formalism. 

The paper is organized as follows: Section two is concerned with the dynamics of the 
FRW cosmology in modified non local gravity. In section three we study the stability of 
the autonomous system of differential equations. We consider a linear perturbation for 
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our model and investigate the evolution of the cosmological perturbations and identify the 
attractor solution for the perturbations during the tracking regime. We also examine the 
dynamics of the universe by stability analysis. Finally, we summaries our paper in section 
four . 

2. THE MODEL 

We start with the action of the non-local gravity as a simple modified garvity given by 
3o|, 

S = l d'x^S^^Ril + fiD-'R))^ (1) 

where Mp is Plank mass, / is some function and □ is d'Almbertian for scalar field. Generally 
speaking, such non-local effective action, derived from string theory, may be induced by 
quantum effects. A Bi- scalar reformation of non-local action can be presented by introducing 
two scalar fields (p and ijj, where changes the above action to a local from: 

S = I d'x^i^^[R{l + f{<P)) + ij{D<p-R)]^, (2) 

where ip, at this stage, plays role of a lagrange multiplier. One might further rewrite the 
above action as 

S = I d'xV^ [R{1 + /(</.) - ^) - d,^pd''^]y (3) 

which now is equivalent to a local model with two extra degrees of freedom. By the variation 
over ip, we obtain = R oi (p = n~^R. 

Now in a FRW cosmological model with invariance of the action under changing fields 
and vanishing variations at the boundary, the equations of motion for only time dependent 
scalar fields, (p and tp, become 

4> + 3H(p + R = 0, (4) 
ij + mi)-Rf = 0, (5) 

where the scalar curvature R is R = 12H^ + 6H, H is Hubble parameter and prime means 
derivative with respect to the scalar field. Variation of action ([3]) with respect to the metric 
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tensor (7^^ gives, 
1 



= 2^?M- {i?(l + / - ^) - dp^pdP<P} - R^,{1 + f-^) + -id^^d,(j) + d^<pd,^lj) 



The 00 and ii components of the equation ([6]) are 



(6) 







1 : : 



d 



= {2H + 3H'){1 + / - ^) + ^0 + ( + 2H-){f - ^) 



2"" ' 'df^ 
Equations and can be rewritten as 

1^'"^ 3Hif'<P-i;) 



dt' 



(7) 
(8) 



3if2 = ^ 
2H + 3H^ = - 



(1 + /-^) ' 
lV'0+(^ + 2iJ|)(/-^) 



(9) 
(10) 



(1 + /-^^) 

Comparison with the standard Friedman equations, the right hand side of the equations ([9 
and (fTOjl can be treated as the effective energy density and pressure: 



M2 

p 

Peff 



(1 + /-^) ' 



M2 



(11) 

(12) 



(1 + /-^) 

Using Eqs. (jlj) and and doing some algebraic calculation we can read the effective energy 
density and pressure from the above as, 

'1 



Peff 



M2 
p 



(13) 



Peff 



M2 

p 



1 + / - ^ - 6f 2 



1 . . . . /' QH{f^-^) 



1 + 



Now by using Eqs. ( IT3l) and ( |T^ the conservation equation can be obtained as, 

Peff + 3H Peff (l+Ueff) =0, 



(14) 



(15) 



where cug// is the EoS parameter of the model. 
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3. PERTURBATION AND STABILITY 



In this section, we study the structure of the dynamical system via phase plane analysis, 
by introducing the following four dimensionless variables, 

x = -f{(f)), y = jj^ ^^6^' ^^^^ 

From Friedmann Eq.([9]) and assuming /(0) = foe^'^, one finds the Friedmann constraint 
equation in terms of the new dynamical variables as 

l = x + k + zy + axy + Qz. (17) 

Then using equations ©-([S]), fl9|)-f lT0|) . and f fTSj) the evolution equations of these variables 
become, 

J doc tic / -* 

k' = -— = -=Qz. 21 
a ma H 

Also, we can obtain effective EoS parameter and deceleration parameter as g = — 1 — 
cUeff = —1 — 1"^! where in terms of the new dynamical variables we have, 

H 4ax{6 + y) +Az — a'^xy'^ + 6zy 

IP^ 2{l-x{l-6a)-k) ^ ' 

For an expanding universe with a scale factor a{t) given by a oc one also can find 

p_ 2{l-x{l-6a)-k) 

Aax{6 + y) + Az — a'^xy'^ + Qzy' 

in terms of the new dynamical variables. We will restrict our discussion on the existence 
and stability of critical points in an expanding universes with H > 0. Critical (fixed) points 
correspond to points where x' = 0, y' = 0, z' = and k' = 0. 

In order to study the stability of the critical points, using the Friedman constraint equa- 
tion (fT7|) we first reduce Eqs. (ITH]) - (PTl) to three independent equations for x, y and z. Sub- 
stituting linear perturbations x x + 6x, y ^ y + 6y and z ^ z + 6z about the critical 
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points into the above three independent equations to first-order in the perturbations, gives 
the evolution equations of the linear perturbations, which yield three eigenvalues Aj. Stabil- 
ity requires the real part of all eigenvalues to be negative. The linearization of the system 
about these fixed points yields three eigenvalues 

Ai--y, A2---, A3---. 

After solving the new set of equations we find that there is only one critical point for our 
system with the physical property presented in table 1. As stated the critical point is a 
stable point in our model for a = — 



Label 


X 


y 


z 




a 


p 


Stability 


S 


1 








1 

3 


1 


1 

2 


stable 



TABLE L The properties of the critical point. 



The stable critical point, S, represents a scalar field-dominated solution in the radiation 
dominated evolutive universe with a oc t^. In Fig. 1, the attraction of trajectories to the 
critical point in the three dimensional phase plane is shown for the given initial conditions. 
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Fig. 1: The attractor property of the dynamical system in the phase plane. 
Initial values are: (blue): x(0) = 1.0099, y{0) = 0.002, z{0) = -0.004. 
(black): x(0) = 1.0099, y{0) = 0.002, ^(0) = -0.003. 
(green): x(0) = 1.0099, y{0) = 0.002, z{0) = -0.002. 

Fig. 2 shows the projection of the three dimensional phase space onto the two dimen- 
sional y = 0.002 phase space in terms of the dynamical variables x and z. 
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Fig. 2: The phase plane of variation of z respect of a;, asymptotically y = 0.002 
stable equilibrium sink. Initial values are: (blue): x(0) = 1.0099, 2(0) = —0.004. 
(black): x(0) = 1.0099, 2(0) = -0.003. (green): x(0) = 1.0099, z{0) = -0.002. 

Fig. 3 shows the projection of the three dimensional phase space onto the two dimen- 
sional X = 1.01 phase space in terms of the dynamical variables y and z. 
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Fig. 3: The phase plane of variation of z respect of y, asymptotically x = 1.01 
stable equilibrium sink. Initial values are: (blue): y{0) = 0.002, z{0) = —0.004. 
(black): y{0) = 0.002, z{0) = -0.003. (green): y{0) = 0.002, z{0) = -0.002. 

As can be seen, in Figs. 2 and 3 all the trajectories for different values of initial conditions 
asymptotically approaches the projected stable critical point in their plane. It means that 
in a dynamical system sense, this critical point is an attractor and the system is drawn to 
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this point in a self-organizing manner. 

At this stage we study the cosmological evolution of EoS parameter, Uefj in terms of 
the new dynamical variables. As can be seen in Fig. 4, for the stable critical point, by 
perturbation, the trajectory of the EoS parameter with the given initial conditions near 
the critical point values, starts from the state in the past with u>eff = 1/3, undergoes 
a quintessence dominated era with —1 < Wg// < sometimes in that past and future, 
passes the current negative value Ueff — —0.14 and finally approaches the same state with 
ujeff = 1/3 in the future, which correspond to the stable critical point (blue curve). In this 
Figure we also show, for completeness, that if the the initial conditions are far enough from 
the critical point values, by perturbation, the universe starts from the state in the past with 
Ueff = 1/3 and undergoes phantom crossing twice in the past and future and never come 
pack to the same state in the future ( the red curve). 




-1 -0.5 0.5 1 

Ina 



Fig. 4: The graph of Ueff respect of ln(a). Initial values are: 
(red): x(0) = -2, y{0) = -2.5, ^(0) = 0.3, 
(blue): x(0) = 1.01, y{0) = -0.04, z{0) = -0.02, 

4. SUMMARY AND REMARKS 

In this paper, we consider a local scalar-tensor formulation of non-local gravity as a simple 
modified model characterized by two scalar fields </> and and function /(</>) which can be 
viewed as scalar potential in the model. 
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We have made a perturbative analysis of the cosmological model, to investigate the 
stability of the model. The perturbative equations of motion are solved numerically, and we 
found that the system is stable for only one critical point under scalar perturbation. The 
three dimensional phase space of the model gives the corresponding conditions for tracking 
attractor. We have also shown the projected two dimensional phase space in x = constant 
and y = constant planes. The stability analysis shows that, for the only stable critical point, 
the phantom crossing never occurs. After perturbation, the universe approaches a stable 
attractor in the future with Ue/f = 1/3 (radiation dominated era) and p = 1/2 (decelerating 
expansion era). The universe may undergoes phantom crossing if the chosen initial condition 
are far enough from the critical point values, such that the universe never come back to the 
initially stable state. 
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